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Abstract 

We derive transformation formulas for the generalized polarization tensors under rigid 
motions and scaling in three dimensions, and use them to construct an infinite number of 
invariants under those transformations. These invariants can be used as shape descriptors for 
dictionary matching. 
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1 Introduction 

Shape of a domain can be represented in terms of various physical and geometric quantities such 
as eigenvalues, capacity and moments. The generalized polarization tensor (GPT) is one of them. 
GPTs are an (infinite) sequence of tensors associated with inclusions (domains) and they appear 
naturally in the far field expansion of the perturbation of the electrical field in the presence of 
the inclusion. They are geometric quantities in the sense that the full set of GPTs completely 
determines the shape of the inclusion as proved in [5]. This suggests that GPTs has richer infor- 
mation on the shape than eigenvalues since the full set of eigenvalues does not determine the shape 
uniquely [TO]- Moreover, recent studies [31 IB] show that we can use a first few terms of GPTs to 
recover a good approximation of actual shape of the inclusion. Even the topology (the number of 
components) can be recovered. GPTs have been used for imaging diametrically small inclusions 
and computation of effective properties of dilute composites. We refer to [B] and references therein 
for these applications. It is worth mentioning that the notion of GPTs has been used not only for 
imaging but also for invisibility cloaking [7]. 

When the domain is transformed by a rigid motion or a dilation, the corresponding GPTs 
change according certain rules and it is possible to construct as combinations of GPTs invariants 
under these transformations. This property makes GPTs suitable for the dictionary matching 
problem. The dictionary matching problem is to identify the object in the dictionary when the 
target object is identical to one of the objects in the dictionary up to shifting, rotation, and scaling. 
The standard method of dictionary matching is to construct invariants, called shape descriptors, 
under rigid motions and scaling, and to compare those invariants, and a common way to construct 
such invariants uses the moments In the recent paper [1] new invariants are constructed 

using GPTs in two dimensions: it is shown that a first few terms of GPTs and hence invariants of 
the target can be computed using the measurement of the multi-static response matrix, and then 
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the dictionary matching technique is applied for target identification. Viability of the method is 
demonstrated by numerical experiments. 

It is the purpose of this paper to extend results of [T] to construct invariants using GPTs in three 
dimensions. In fact, using contracted GPTs (CGPT), which are harmonic combinations of GPTs, 
we are able to construct an infinite number of shape descriptors which are invariant under rigid 
motions and scaling. Since the method of above mentioned paper uses the complex structure of the 
two dimensional space, it can not be applied to three dimensions. Using transformation formulas 
of spherical harmonics under rigid motions, which can be found in |13j for example, we are able 
to derive transformation formulas obeyed by CGPTs. We then use these formulas and method of 
registration (see [Tl]) to construct invariants which can be used for target shape description and 
position and orientation tracking [TJ [5] . 



2 Neumann-Poincare operator and CGPT 



For a given bounded domain D in R 3 with the Lipschitz boundary, the Neumann-Poincare operator 
JC d for a density function <fi g L 2 (dD) and its L 2 -adjoint K* D are defined in the principal value by 



K, D [4>](x) = 

47T Jg 



J_ f (y~x,v(y)) 

47r JdD \x - y\ 3 

x - y,v(x)) 



dD k-s/l £ 



4>{y)d(j{y), xedD : 
(p(y)da(y), x € dD . 



Here (■, ■) denotes the scalar product and v(x) the unit outward normal vector along the boundary 
at x. Let A be a real number such that |A| > 1/2. The generalized polarization tensor (GPT) M a p 
for multi-indices a — (ai, a.%, a^) and ft = (/3i, 02, Ps) associated with A and D is defined by 



M afi (X,D)= f /(A7-/C^)- X [i/- Wy a ]da. 

JdD 



(2.1) 



Here y a = y^y^y^ 3 - Throughout this paper A is fixed, so we use M a p(D) for M a p{\, D). 

The contracted GPTs (CGPT) arc harmonic combinations of GPTs. To be more precise, let 
Y™, — n < m < n, be the (complex) spherical harmonic of homogeneous degree n and order m, 
i.e., 



Yl 



1<P) = (-!)' 



2n + 1 (n — m) ! 



1/2 



e imtp V™{cos6), ~n<m<n, 



An (n + m)\ 

where V™ are the associated Legendre polynomials of degree n and order m. If 



then CGPT M mnkl is defined by 



In other words, we have 



ilk 



\ct\—n 



a™ n a k p l M a p, m,n,k,l = 1,2,- ■ 



M nmlk = / r^y ; fe (^,^)(A/ - K* D y 

JdD 



-r-Y™{e yiV y) 



dD 



(y)do{y), 



(2.2) 



(2.3) 



where y = r y (cos <p y sin 6 y , sin ip y sin 9 y , cos 9 y ) . 

We now show that the operator JC* D is invariant under the shift, scaling, and rotation of the 
domain £>, which will be a crucial fact in study of invariancc properties of CGPTs in following 
sections. Let z be a point in 1R 3 , s a positive number, and R a 3 x 3 orthogonal matrix, and define 
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• Shift: T Z D = D z = {x z = x + z \ x £ D} ; 

• Scaling: sD = {sx \ x € D} ; 

• Rotation: R(D) = D R = {x R = Rx \ x € D} . 
For a function (f> defined on dD, define 

4> z {y z ) = 4>{y), s ( S y) = 0(y), <f> R (y R ) = 4>{y) 
Then one can see easily that the following invariance holds: 



K* Dz W\{x z )=K* D [<j>]{x) (2.4) 
K* sD [^](sx)=]C* D mx) (2.5) 

K*D R [<t> R ](xR)=lC* D [<t>}(x)- (2-6) 



In fact, for example, we have by the simple change of variables y = y/s 

ir* ^ 1 f (sx-y,v(sx)) , 

4tt J d(sD) \sx - y\ 6 



d{sD) \sx-y\ s 

(f>(y)s 2 d<r(y) = fC* D [4>](x) 



1 f s(x- y,v(sx)) ± 



4tt J 9d s 3 |a;-?;| 3 
The other two relation can be seen similarly. 

3 The MSR matrix and the CGPT block matrix 

We will investigate invariance of the block matrix consisting of CGPTs rather than individual 
GPTs M a p. The way to construct the block matrix can be seen most clearly using the multi-static 
response (MSR) matrix. 

Let {x r }^ =1 and {x s }^ =1 be a set of electric potential point detectors and electric point sources. 
Let u s (x) be the solution to the transmission problem in the presence of inclusion D 

' V-(1 + (k-1) X d)Vu s (x) = 6 Xs {x), xeR 2 \dD, 

u s (x)\+ = u s {x)\- , xedD, , 

v{x) ■ {Vu s )\ + = kv{x) ■ (Vua)|- , xedD, { ' 

u s (x) -T s (x) = 0(\x\- 2 ) |x-x s |^oo, 

where the notation 4>\±(x) means the limit lim t _>.o (j){x ± tu(x)), n = (2A + 1)/(2A — 1), and 

1 1 



r s (x) = r(x - x s ) 



47T \x — I, 



Without the inclusion D the solution is r s (x). The MSR matrix V is the matrix of differences of 
electric potentials with and without the conductivity inclusions, and its rs-components are defined 
by 

V rs = u s (x r ) - T s (x r ), l<r,s<N. 
Let Sd be the single layer potential associated with D: 

S D [<f\ (x) := / r(x - y)ct>(y)da(y) , x G R 3 . 

JdD 



It is known that the solution u s to (|3.1[) can be represented as 

u s (x) = T s (x)+S D [^ s }(x) (3.2) 
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where 4> s £ L 2 (dD) solves 



So, the MSR matrix is given by 



(\I-K* D )i<t> a ] = ^ 
ov 



i)D 



v„ = 



\ Y[x r -y){\l-K* u )- x 

JdD 



~dr s 




dv 


dD. 



(y)da(y) 



Let x = r x (cos tp x sin 9 X , sin ip x sin 9 X , cos 6 X ) and y = r y (cos ip y sm.0 y , sinip y sm.Q y , cos 6 y ) in 
spherical coordinates, and suppose that r y < r x . Then it is well known (see [I21II2JEI for example) 
that 

i=o k=-l Tx 
So, by assuming the inclusion D is away from the sources, we have from (|2.3p 

oo l ^ i 



(3.3) 



r / 00 1 r i \ 

^ = / flD (E E ^(•.^w^)^ 
\ l=0 k= _ l 

x (A/-/CB)- 1 



00 n 



00 n 



EEE E 



£=1 fe=— 1 n= 1 m— — n 



^E E ^^(WW^J^r 

n— m— — n s 

(2Z + l)(2n+l) ( ^ ' ^ } 



<7D 



i.v)d<y{y) 



x( f r[Y^9 y ^ y ){\I-K,* D )- l \^Y: n {e y ^ y ) \y) d( x{y)]4+T Y ^ Xs ,V Xs ) 

\JdD Y 0v dD\ / r z 3 

= EEE E f^^ ^^t^^jMn^^^j. (3.4) 

Here we have used the fact that M nm ik = for n = or Z = 0. 
We now introduce matrix M;„ by 

(M ; „) fem := M nm i k , -Z < fc < Z, -n < m < n. (3.5) 

We emphasize that the dimension of M;„ is (21 + 1) x (2n + 1). We also define 1 x (21 + 1) and 
1 x (2n + 1) matrices (vectors) Y r j and Y sn by 



(Y sn )m 



(2Z + lHt 1 : 



-l<k<l. 



-n < rn < n . 



(3.6) 
(3.7) 



(2U+1K+ 1 ' 

Then (|3.4p yields, after truncating terms corresponding to I > A' and n > K for some integer if 

K 



F sr = £ Y rl M ln Y* n , 



(3.8) 



2,n=l 



where * denotes the Hermitian transpose (^4* = (AY, t for transpose). 
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Let N be the number of (coincident) source transmitters and receivers and K be the truncation 
order as before. We further define 



M 



Mn M12 
M 2 i M 22 



Mijc 
M 2K 



Mjfi M K2 ■■■ M 



KK 



Y11 Y12 
Y21 Y22 



Y lK 
Y 2K 



Y N1 Y 



N2 



Y 



NK 



(3.9) 



Then M and Y are (K 2 + 2K) x (K 2 + 2K) and N x (K 2 + 2K) matrices, respectively, and the 
MSR matrix V can be written as 

V = YMY* (3.10) 

after truncation. We call M;„ and M a CGPT matrix and the CGPT block matrix of order K, 
respectively. 

Proposition 3.1 The CGPTs matrix M is Hermitian, i.e., M = M*. Furthermore, the matrices 
M„„ are invertible for all n > 1. 

Proof. Define the coefficients a„ and a™ n so that 

\P\=l \a\=n 

Then from (|2.2p and the symmetry property of Af Qj( g on the coefficients of harmonic polynomials 
13 Theorem 4.10], we have 



M nmlk = J2 a k l (C n M af> = af< n Mp a = M lknm , 

\a\=n,\fi\=l \a\=n,\0\=l 

which yields M ; „ = M* l . So M is Hermitian. 

To prove the invertibility of M„„, it suffices to show v*~M. nn v 7^ for any v € C 2n+1 , b ^ 0. 
By definition of the CGPT matrix, we have 

2n+l 2n+l 
V*M nn V = X! X! V *k M nrnnkVrn 
k=l m=l 



r 2n+l r „ 2)1+1 

/ E Vkr?Y*(0 v , y y ){\I VnT 1 — v m r^Y™(6 y ,Vy) 

j 9d fe= i yov m=i 



dD. 



(y)dcr(y) . 



Note that 2feS v kfy Y„ (Qy: fy) i s a harmonic polynomial, and we introduce the coefficients j a 



so that 



2n+l 



E W = E v k r^{9 v ^ y ). 

\a\—n 



k=l 



And write 7a = 7 Q + Z7 Q , with 7* and 7* the real and imaginary part, respectively. Remark that 
both 7* and 7^ are coefficients of some harmonic polynomials. Then using the fact that u*M„ n « 
is a real number, we get 

v*M nn v= 1* M *P1P= E Ta M «^+7 a M^ 7 |, 



M = 

101= 



l«l= 

101 = 



while this quantity is strictly positive if A > 1/2, and strictly negative if A < 1/2 by the positivity 
of M a p [6j Theorem 4.11]. This completes the proof. □ 
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4 Transformation formulas for the CGPT matrix 



In this section, we derive transformation formulas of the CGPT matrix M;„ defined in (|3.5|) and 
the CGPT block matrix under rigid motions and dilation. These formulas play an essential role in 
finding the invariant under these transformations in the sequel. 

4.1 Scaling 

We consider first the scaling of M/„. Let us denote 



4>D,nm{y) = (A - K* D )~ 



(y)- 



dDi 

Using the change of variables y s = sy (r Va = sr yi 9 Vs = 9 y , ip ys = tp y ), we obtain 
M nmlk (D) = / r l Y k l (9 y ,ip y )(j) D ^ m (y)da(y) 

JdD 

s ~ lr y Y k [ (9 ys ,<p ys )<f) D . nm (-y s )s~ 2 do-(y s ) 

d(sD) s 



r y y k (0y a , f ys WD,nm(ys)dcj(y s ) . 

d(sD) 



Since 



it follows from (12.51) that 



(XI - KdWd^Vs) = (XI - fC* D )[<f> D , nm ](y) 
= (u(y),V(r-Y-(e y ^ y )) 
= S - n+1 (v(y s ),V(rlY r l(9 ys ^ ys )) , 

and hence 

<t> s D, nm {ys) = s~ n+1 <t> S D,nm(y s ) ■ 

Thus, we have 

M nmlk (D) = / rl a Y k \e ys ,ip ys )4> sD>nm (y a )da(y s ) = s -«+ n +V M nmlk (sD) . 

JdsD 

Lemma 4.1 (Scaling) For any positive integers l,n and scaling parameter s > 0, the following 
holds: 

M ln (sD) = s l+n+1 M ln (D) . (4.1) 

4.2 Shifting 

To deal with a shifting of M;„ , we need the translation of the regular spherical harmonics r n Y™ (9,ip). 
For y = (r, 9, (p), z = (r z , 9 Z , ip z ), and y z = y + z = (r' , 9' , ip') , we have the following expression of 
the translation of the regular spherical harmonic, 

(n,m) 

r' n Y™(6', ^ = ]T C^ nm rr v Y™-S(0 z ,t Pz )r»Y l ?(6, <p) , (4.2) 



where 



fj,nm 



47r(2n + l)(n - m)!(n + to)! 



(2n -2v + l)(2v + l)(n -v-m + fi)\(n -u + m- n)\{y - n)\(u + 



1/2 
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Here we use special summation notation: 

(n,m) n ^l— min(i/, — f+n-j-m) 

E = E E ■ 

(i>,aO i/— /i— max( — n+m) 

We refer to [13] for a proof of formula (|4.2|) . Using (|2.4|) and (|4.2|) . we obtain 
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M nmlk {D z ) = / ^(^^(AZ-K?,)- 1 
Jan 



dD. 



(y)da(y) 



r il,k) 



x (\I-K,* D )- 1 



(n,m) 



dv 



E C^ nrn r n z - u Y™-/(6 z ,tp z yY»(6, <p) 



WD. 



{l,k) (n,m) 



(4.3) 



For a pair of integers i and /, define (21 + 1) x (2i + 1) matrix G;j = Gu(z) by 

Ci 2 ihr l ~ % Y^ (9 z ,ip z ) if max(-i, i - / + k) < j < mm(i, -i + l + k), 



(Gli)kj — Gijlk 







otherwise 



(4.4) 

for —l<k<l. Then, from the equation (|4.3[) and taking the fact that M;„ = if n . — or I = 
into account, we immediately have the following property. 

Lemma 4.2 (Shifting) For any positive integer l,n and the shifting parameter z, the following 
holds: 

l n 

MhfD^^G^jM^fl)^^. (4.5) 
i=i v= \ 

To gain better understanding of the matrix Gu, we compute G21. Indeed it will play a role in 
Section [SJ G21 is given by 



Go 



For k = —2, —1, 0, 1, 2 (I = 2 and i = 1), the conditions max(— i, i — I + k) < j < min(i, — i + I + k) 
can be written as follows: 



Gi,- 


1,2,-2 


Gi.0,2,-2 


"1,1,2,-2 




1,2,-1 


Gl, 0,2,-1 


"1,1,2,-1 


Gi, 


-1,2,0 


G\, 0,2,0 


"1,1,2,0 


Gi, 


-1,2,1 


"1,0,2,1 


"1,1,2,1 


Gi, 


-1,2,2 


"1,0,2,2 


"1,1,2,2 



k = 


-2 


max(- 


1, 


-1 - 


hfc) 


<3 


< 


min(l, 1 - 


-k) 


=> 3 


= -1 


k = 


-1 


max(- 


1, 


-1 - 


hfe) 


<3 


< 


min(l, 1 - 


-k) = 


>3 = 


-1,0 


k 


= 


max(- 


1, 


-1 - 


hfc) 


<3 


< 


min(l, 1 - 


-k) = 


>3 = 


-1,0,1 


k 


= 1 


max(- 


1, 


-1 - 


Vk) 


<3 


< 


min(l, 1 - 


-k) = 


>i = 


0,1 


k 


= 2 


max(- 


1, 


-1 - 


Vh) 


<3 


< 


min(l, 1 + k) = 


>3 = 


1. 



r 



Then using (|4.4[) and the definition of the spherical harmonics, we have 





~Gi,- 


-1,2,-2 











G h - 


-1,2,-1 Gi,o,2,-l 









G21 — 


Gx, 


-1,2,0 Gi : o,2,0 


Gi, 1,2,0 










Gi. 0,2,1 


Gi, 1,2,1 













Gl,l,2,2_ 








"Gi,- 


i,2,-2r z Yf 1 (O z ,ip z 


) 









Ci, 


-i,2,-ir z Y?(6 z , Vz ) 


Gi.o^.-i^Ff 1 


(Ox,<Px 




Ci 


-i,2flr z Y^(9 z ,tp z ) 


Gi.o,2.o^n°(^ 


] z,Vz) 









Cx A 2,ir z Y£{t 



















Cl, 1.2,0^^1 1 (0z, ( Pz) 

Ci,i,2,ir z Y?(6 z ,<p z ) 
Ci,i X2 r z Yl{e z ^ z ) 



-VZr z sin z e- itp * 
\fhr z cos 8 Z 



%r z siii^e 1 ^ 





r, sin 9 z e~ 



^■r z cost 










^r z sin 6 z e~ l 

\fhr z cos Z 
sin O^f* 



(4.6) 



On the other hand, from the equation (|4.4[) . one can easily see that the G n jnk is non-zero only 
when fc = j, and 

G nfcnfc = Cnfcnfcltf = ^Y °(9 z ,<p z ) = 1 . 

Thus, G„„ is the identity (2n + 1) x (2n + 1) matrix. So, (|4.5[) yields, for instance, 



M 21 {D Z ) = G a i(*)Mu(D) + M 2 i(£>) 



(4.7) 



4.3 Rotation 



Rotations in three dimensions may be described in many different ways. Among them we use 
the Euler angles which can be conveniently used to represent the rotation formula for spherical 
harmonics. The rotation R is given by 



R 



cos 7 — sin 7 
sin 7 cos 7 
1 



cos /? — sin (3 

1 
sin /3 cos f3 



cos a —sin a 
sin a cos a 
1 



That is, we rotate by an angle a about z-axis, and by an angle (3 about the new y axis, and 
finally by an angle 7 about the new z-axis. Let R be a rotation matrix. Since the homogeneous 
polynomials and Laplace operator are invariant under rotation, Y™ (i?£) is also spherical harmonic 
of degree n, moreover Y™(R£) can be written as follows: 



where 
Here. 



XT (0 , 

m',m im' y jm' ,m / a\ ima 

Pn ~ e a n \P) e ! 

l {/3) = [(n + m')l(n - m')\(n + m)!(n - m)!] 1/2 

^_ l) m '+ m + s (cos ^) 2 (™~ m '(sin ^) 2s ~ m + m ' 



(4.8) 
(4.9) 



(n + m — s)!s!(m' — m + s)!(n — m' — s)! 
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where the sum is over values s such that the factorials are nonnegative: 

max(0, m — m ) < s < min(n — m! , n + m) . 

A proof of (|4.8[) can be found in [15] . 

We want to find the transformation formula M nm n~, i.e., M nm ik(Dfi). We have 



M nmlk (D R )= f r% k {RO{XI-JC* DR ) 

JdD„ 



f)D r , 



(Ry)da(Ry) 



[ r l J2 pf' fe ^ fe '(e)(A/-/CB)- 1 [|-r" Pn' m Y™'(® 

JdD i, — , V au 7 — 

k'— — l 

qFM ; „( q ™) 4 , 



m — — n 



(4.10) 



where 



n,m n,m\ 
n 5 * * * 5 rn ) 



Let for each positive integer n 



Qn — Qn{R) 



~—n,—n 



„— n, — n+1 „— n+1,— n+1 
"n rn 



-n+l,n 



rn 
n n, — n+1 



(4.11) 



r n In rn 

The rotation matrix Q„ is called Wigncr _D-matrix and known to be is unitary (see [13]). 
Lemma 4.3 (Rotation) For a unitary matrix R the following relation holds: 



M ln (D R ) = Qi{R)Mi n {D)Q(R) t 
4.4 CGPT block matrices and its properties 



(4.12) 



In addition to the CGPT block matrix, we define two block matrices: Let K be the truncation 
order, s be a scaling parameter, z shifting factor, and R rotation, and define 



G(R) := 



Gn 
G21 G22 



Gki G 



A' 2 



G 



A' A' 



, Q(s,R) 



sQi 

.s 2 Q 2 











A 



(4.13) 



where G;„ = G;„(z) and Q„ = Q n (R) are defined by (|4.4[) and (|4.11[l . respectively. Then we have 
the following theorem 

Theorem 4.4 Let T z be the shift by z, T s scaling by s, and R a unitary matrix. Then the CGPT 
block matrix M satisfies 



M(T Z T S R(D)) = sG(z)Q(s, R)M(D)Q(s, i?) 4 G(z)' 
Proof. We have from (|4.5[) that 

M(T Z T S R{D)) = GM{T s R(D))G t . 



(4.14) 
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We also have from (|4~T|) and (|4.12[) that 

s 2 M u (fi(i})) 



M(T S R{D)) = s 



So we obtain P~T4l . 



a K + l M K1 (R(D)) 
s 2 Q7M u (D)Q* 



S 1+A 'M 1 ^(i?(D)y 



^Mjfjf^D)) 
s 1+if QTM 1K (£)Q 



2f 



□ 



5 Transform invariant shape descriptors 

In this section, we construct the invariants using CGPT matrices under shifting, scaling, and 
rotation. 

Let B be a reference domain and D be the one obtained by rotating B by R, scaling by s, and 
shifting by z, i.e., 

D = T Z T S R{B). (5.1) 
Since Q„ is unitary, we have Q„ 1 = Q^. One can easily see that Mix an d Qi commute, i.e., 

M11Q1 = Q1M11. (5.2) 
Since G nn = I 2 n+i ((2w + 1) x (2n + 1) identity matrix), we have 
Mu(X>) = s 3 Q7M 11 (B)Qi = s 3 M n (B) , 



M 21 (D) = s J G 21 Q 1 M 11 (B)Qi + s 4 Q 2 M 21 (B)Q t 1 = G 2 iM u (Z>) + s 4 Q 2 M 21 (B)Qi 



Define 



Then we have 



U D = M 21 (D)M 11 (D)- 1 and U B = M 21 (B)M n (B)' 1 . 



U D = G21 + ,sQ 2 U B Q t 1 . 



(5.3) 



5.1 Invariance by registration 



We first recall that the matrices G;„ and Q„ are determined by the shift factor z and rotation R, 
respectively. Moreover, in view of (|4.6p . we may write G 2 i in terms of rectangular coordinates as 



G 21 (z) = VE 



(zi + z 2 i) 
Z3 

§0*1 ~ z 2 i) 









+ z 2 i) 



3 z 3 

\{z\ - z 2 i) 







\{z\ + z 2 i) 



Z3 

z\ - z 2 i 



(5.4) 



We present here a method of registration to construct invariants. This method is based on a 
linear mapping u : A/5 x 3(C) — s- C 3 (Ms X 3(C) is the collection of 5 x 3 complex matrices) which 
satisfies 

u(G 2 i(z)) = z (5.5) 
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and 



u(Q 2 (i?)UQi(i?) t ) = Ru(XJ) for any U G il/ 5x3 (C) and any rotation R. 
Such a linear mapping docs exist and an example is 



1 



" (U) " 7t 



^ Ul1 + 102 U31 _ I572 U22 + T^f U42 - 10: 

3 TT I TT _L 3 TT L 3 TT L TT__ I 3 



l (lO U ll + T^ U 31 + W2 U 22 + I57f U 42 + Y^Tf U 33 + YqUss) 

/T 



^u 21 + fu 32 + iu 43 



(5.6) 



(5.7) 



One can easily check that this linear transformation satisfies (|5.5p using (|5.4j) . We check that it 
satisfies (|5.6I) through symbolic computations using Matlab. Since the computation is lengthy, we 
omit the detail here. 

In the following we use the shorthand ud = w(Ud). Applying u to (|5.3[) . we have 



ud = z + sRub, 



(5.8) 



somc can be used as a registration point (see [H]). We emphasize that since = M 2 i(-D)Mn(Z)) _1 , 
Ud can be computed using CGPTs of D. 
The first invariant we introduce is 



J ln (D) := M ln (T_ UD D) = Gi{-u D )M ln (D)G l (-u D ) t l,n= 1,2,., 
Proposition 5.1 For all indices l,n the quantity J\ n is invariant by translation: 

Jin(T z D) = J ln (D) for any z£l 3 . 
Proof. Let D z — T Z D. By definition of Ji n , and using formula (|4.5[) . we have 

Jin{T z D) = M ln {T_ UD D z ) = M ln (T z _ UB D) 



(5.9) 
(5.10) 



= Gj(z - «„JM,„(D)G„(z - u D J*. 
The relation f|5 . 8[) says that 11^. = z + ud, and hence we have 



Jin{T z D) = G l (-u D )M{D)G n (-u D ) t = J ln (D). 



This completes the proof. 

In the following we denote D = D — ud, then 



□ 



Lemma 5.2 Let D = T— UD D. For any indices l,n and scaling parameter s > and rotation R, 
we have 



J ln (sR(D)) = Ji n (sR(D)) = M ln (sR(D)) 



J+n+l 



M ln (R(D)). 



In particular, we have 
and 



J ln {sD) = s l+n+L Ji n (D) 



J ln (R(D)) = QiWJtoMQnffl. 
Proof. Since J\ n is translation invariant, we have for any z G M 3 : 

J ln {sR(D)) = J ln (T sRz sR(D)) = J ln {sR(T z D)). 
Then by taking z = -tic, we obtain the first identity in (|5.11[) . Moreover, we have from ([5 

= — ud + u D = 0. 



(5.11) 
(5.12) 
(5.13) 
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So we have u(sRD) = sRu^ = , which implies the second identity in (|5.11[) . The third one is 
(|4~Tj) . By taking R = I, we have (|5T2|) . (|5T3)) follows from the definition of J ln and (|4T2|) . □ 
In particular, it can be seen from (|5.11[) that all J nn are invertiblc. So we define the second 
invariant: 

Sln(D) := JuniD^Jnl^JuiD)^^^). (5.14) 

It is worth emphasizing that Si n (D) is a square matrix of dimension (2n +1). 

Proposition 5.3 For any indices l,n, the quantity Si n {D) is translation and scaling invariant: 

Sin{T z sD) = Si n {D) for any z £ R 3 , and s > 0. (5.15) 
Moreover, for any rotation R: 



S ln (R(D)) = Q n (R)S ln (D)Q n (RY . (5.16) 

Proof. Since J^n's are translation invariant, so are Si n 's. Scaling invariance of Si n follows from 
()5.12[) . Formula (|5.16[) can be seen using (|5.13[) as follows: 

S ln (R(D)) = J nn {R{D))- 1 J nl {R{D))MR{D))- 1 J ln {R{D)) 
= Q^Jnn(D)- 1 J nl (D)MD)- 1 J ln (D)Q t n 

This completes the proof. □ 
Since the Frobenius norm of a matrix remains unchanged after the multiplication by a unitary 
matrix, so finally we define the shift descriptor 

X ln (D) = ||5, n (D)||jp (5.17) 

which is clearly invariant by any rotation, scaling and translation. 



6 Concluding remarks 

In this paper we have constructed new shape descriptors in three dimensions which arc invariant 
under translation, rotation, and scaling. These shape descriptors can be used to efficiently identify 
a target using a dictionary of precomputcd CGPTs data. They can be also used for tracking the 
position and orientation of a mobile three-dimensional target from multistatic measurements. 
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